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ABSTRACT 



We will show the importance of coexistence of oppositely flowing ^-currents 
for magnetized stars to sustain strong toroidal magnetic fields within the stars by an- 
alyzing stationary states of magnetized stars with surface currents which flow in the 
opposite direction with respect to the bulk currents within the stars. We have imposed 
boundary conditions for currents and toroidal magnetic fields to vanish outside the 
stars. It is important to note that these boundary conditions set an upper limit for 
the total current within the stars. This upper limit for the total current results in the 
presence of an upper limit for the magnitude of the energy for the toroidal magnetic 
fields of the stars. If the stars could have the toroidal surface currents which flow in 
the opposite directions to the internal toroidal currents, the positively flowing inter- 
nal toroidal currents can become stronger than the upper limit value of the current for 
configurations without surface toroidal currents. Thus the energies for the toroidal 
magnetic fields can become much larger than those for the magnetized stars with- 
out surface toroidal currents. We have also analyzed the same phenomena appearing 
in spherical incompressible stars for dipole-like magnetic fields with or without sur- 
face toroidal currents by employing the zero-flux-boundary method. We have applied 
those configurations with surface toroidal currents to magnetars and discussed their 
flares through which magnetic helicities could arise outside the stellar surfaces. 
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1 INTRODUCTION 

The origin of the magnetic fields inside the star is not still un- 
derstood well. According to the recent progresses of the the- 
oretical researches, there are two possibilities for the mech- 
anism to generate and sustain magnetic fields, i.e. 1) the dy- 
namo theory and 2) the fossil field theory (see the review by 
IMossI I1994). Stars which have convective regions regenerate 
and sustain the magnetic fields by the dynamo mechanism. On 
the other hand, A, B, O star s have a convect ive core where 
dynamo action develops (e.g. lBrun et alj|2005l ). However, the 
timescale for the dynamo -generated field in this core to reac h 
the surface is too long l lCharbonneau & MacGregorl 2001). 
Therefore the magnetic field observed at the stellar surface 
is supposed to be fossil field. Recent observations show the 
obs erved magnetic fie lds' properties confirm this hypothesis 
(see lWade et alj201 ll) . Consequently, it has long been consid- 
ered that magnetic fields of such stars would come from fossil 
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magnetic fields. According to the fossil field theory, the mag- 
netic fields of such stars would be originated from magnetized 
interstellar media. If the magnetic fluxes would be conserved 
during star formation processes, the magnetic fields could be 
concentrated to smaller regions of the stars by the gravita- 
tional contraction and result in strong magnetic fields inside 
the stars. Since the electric conductivities of the stars are very 
large, those magnetic fields do not diffuse during their forma- 
tion stages. As for neutron stars, the origin of their strong mag- 
netic fields is much more uncertain. If the fossil field theory 
could be applied, the strong magnetic fields could be reached 
by the same mechanism as those mentioned above for stars 
without external convective regions. On the other hand, if we 
could adopt the dynamo theory, the strong magnetic fields 
would be formed by the dynamo due to the rapid differen- 
tial rotations in the convective regions of the proto-neutron 
stars (Duncan & Thompson 1992). The differential rotations 
could wind up the initial poloidal magnetic fields to produce 
strong toroidal magnetic fields before the crusts would crystal- 
lize. In either case, the magnetic fields of neutron stars would 
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be present even at birth and survive in much longer timescales 
than the Alfven timescale (= ^/4npr 3 /B ~ 100 s for typi- 
cal neutron stars with magnetic fields of order of B — 10 12 G). 
Here r s is the radius of the stellar surface which is, for axisym- 
metric configurations, a function of 6 of the polar coordinates 
(r,9,ip), i.e. r s (0). In order to sustain these kinds of fossil 
magnetic fields for a long time, these magnetic fields must be 
stationary and stable. Thus, in order to understand the mag- 
netic fields originated from the fossil fields, it would be useful 
and important to get stationary and stable configurations of 
magnetized stars. 

As for stability, analytic studies have shown that any con- 
figurations with either p urely poloidal or purely toroidal ma; 
netic fields are unstable (Tayler 1973; Markev & Tavlerll973l 
Stable magnetized stars should have both the poloidal and 
the toroidal magnetic fields. Moreover, the magnitudes of the 
toroida l fields mu st be comparable with those of the poloidal 
fields dTavlenl 19800 . This argum ent has been shown to be the 
case from the recent simulations. Braithwaite & Spruit (2004) 
have shown that an initial random magnetic field in stably 
stratified stellar layers relaxes on the stable twisted-torus mag- 
netic field configuration after several Alfven timescale. Simi- 
lar twisted-torus magnetic field configurations have been also 
obtained in many previous works by numerically exact com- 
putations of the axisymmetric stationary states of magnetized 
starsjJTomimura & Eriguchi 2005; Yoshida & Eriguchi 2006; 
Lande r & Jonesl 120091 ; lLander et al.l 1201 lc iFuiisawa et al.l 

20121) . structure separated Grad-Shafranov (GS) solving 
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method dCiolfi et al. 2009; Glampedakis et al. 2012) or zero- 

flux-boundary method dPr endergast 1956; Io ka & Sasakil 



|2004 iDuez & MathijbOld ; lYoshida et alj|2012h . The stabil- 
ities of these fields, however, have not been clarified yet be- 
cause it is difficult to analyze their stability by linear stabil- 
ity analyses or other mea ns based on the sta tiona ry configu- 
rations . On the other hand, iBraithwaitel (120091) and lDuez et alj 
(2010) have shown that the stability criteria of the magnetized 
stars could be expressed as below: 
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(1) 



where M/\W\ is the ratio of the magnetic energy to the grav- 
itational energy (see Appendix lAU . A4 P /M is the ratio of 
the poloidal magnetic energy to the total magnetic energy 
and a is a certain dimensionless factor of order 10 for main- 
sequence stars and of order 10 3 for neutron stars. The value 
of .M/|W| is about 10~ even for magnetars and is expected 
to be M/\W\ < 10~ 6 for other real stars. Thus the left hand 
side of this inequality could be less than about 10 -3 even if the 
value of a might be ~ 1000. Therefore, this criterion means 
the configurations with the twisted torus magnetic fields are 
stable even if the toroidal magnetic fields are much stronger 
than the poloidal magnetic fields. In contrast, the right hand 
side of this inequality means that the strong poloidal magnetic 
field configurations are unstable (we define that the poloidal 
fields are strong when M p /M > 0.8 in this paper). As 
shown in dynamical simulations mentioned above, configura- 
tions with the strong poloidal magnetic fields are likely to be- 
come unstable within several Alfven timescales. This criteria 
would not be applied to all situations because we might be able 
to consider various kinds of magnetic field configurations as 
the initial states and different choices of the initial conditions 



might influence on the evolutions of the magnetic fields. How- 
ever, it seems to be the case that there is a tendency to become 
more unstable even for the twisted-torus magnetic field con- 
figurations with larger poloidal magnetic field energies. There- 
fore, it would be a natural consequence to consider that there 
would be stable magnetized stars with strong toroidal mag- 
netic fields which satisfy the condition M p /M < 0.8, i.e. 
Mt/M > 0.2 where Mt is the energy of the toroidal field. 

On the other hand, the majority of investigations in 
which stationary states of the magnet i zed stars have been 
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treated (e.g. lYo shida & Eriguch i! l2006t JYo shida et al 
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Lander & JonesP2 009; Cio lfi et aUkOQgJ lLander et al 

Fujisawa et al. 2012) failed to obtain configurations with 
strong toroidal magnetic fields. In these studies, stationary 
states of magnetized stars hav e been pursued either by nu- 
merically exact methods (e.g. iTomimura & Erig uchi 2005) 
or structure separated GS solving methods (e.g. Iciolfi et al.l 
2009). However, they have only found that it was very difficult 
to obtain stationary states of magnetized stars with very strong 
toroidal magnetic fields. In some of their solutions the toroidal 
magnetic fields have been almost as strong as the poloidal 
magnetic fields only in the particular local regions inside the 
stars, but the total energies of the toroidal magnetic fields as a 
whole are much smaller than those of the total poloidal mag- 
netic fields. In other words, the ratios of A4 P /M in their so- 
lutions are much bigger than 0.8. 

By contrast, some studies of magnetized stationary 
configurations by structure separated GS solving method 

(Glampedak is et al.| 120121) by zero-flux-b oundary method 
Duez & M athij l20ld : lYoshida et alj |2012|) have succeeded 
in obtaining the magnetized equilibria with strong toroidal 
magnetic fields by choosing very special boundary condi- 
tions for the pol oidal magnetic fields . The boun dary condi- 
tion ad opted bv lTomimura & Eriguchi! (2005) and lCiolfi et al.l 
(2009) in which they failed to obtain configurations with 
strong toroidal magnetic fields is that the poloidal magnetic 
field lines should continue smoothly through the stellar sur- 
faces into the vacuum region which is considered to be outside 
of the star s. On the other ha nd, the boundary condition em- 
ployed by Glampedaki s et alj J2012I) is different. The poloidal 
magnetic field lines need not continue smoothly at the stellar 
surfaces, because in some of their models it has been allowed 
for the surface currents to exist. By specifying such a boundary 
condition, they have succeeded in finding that magnetized con- 
figurations whose total energies of the toroidal magnetic fields 
become much stronger as the surface currents are increased. 

IDuez & Matiul feoicl) and lYoshida et alj 1201 2h also ob- 
tained the stationary configurations with strong toroidal mag- 
netic fields, but the boundary condition which they adopted is 
of different kind from those mentioned above. Their assump- 
tions are essentially the same as those in the previous works 
JPrendergastl Il956l ; IWoltierl Il959all5 Il960l; lloka & Sasakil 
I2004T ). They imposed the boundary condition that the magnetic 
flux on the stellar surfaces should vanish, so all of poloidal 
field lines are closed and confined inside the stars and no 
poloidal magnetic fields penetrate to the vacuum region out- 
side of stars. In their solutions, the region where the toroidal 
magnetic fields exist inside the star is much larger than that 
of any other models. However, they did not explain the reason 
why the magnetized stars can sustain such configurations with 
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large toroidal magnetic energies under their special boundary 
condition. 

In this paper we will deal with magnetized configurations 
with large amount of the magnetic energies in the toroidal 
fields and present the reason why the magnetized stars can sus- 
tain strong toroidal magnetic fields within the stars. As will be 
shown, we have found that the total currents of the magnetized 
stars are important keys to understand this problem systemat- 
ically and the values of the total currents seem to be deeply 
related to the boundary condition of the magnetic fields. 

It should be noted that for the stationary configurations 
the magnetic fields are governed by the Grad-Shafranov equa- 
tion which is of the elliptic type partial differential equation for 
the magnetic flux function. Therefore, the solutions of the GS 
equation are necessarily strongly depending on the boundary 
condition(s). 

In this paper, w e use both the numerical ly exact 
non-force-free method (ITomimura & Eriguchil 120051) and the 
method in which boundary conditions are applied at finite 
locations from the stellar centre (e . g., [P rendergastj ll956t 
lloka & Sasaki |2004 iDuez & Matrul |201Ch . We will show 
configurations with negative surface currents or with the re- 
gions where the current become negative can sustain the strong 
toroidal magnetic fields inside the star. Here the term 'nega- 
tive' means that the currents flow in the opposite direction to 
the flow direction of the bulk of the interior currents. 



2 FORMULATION 

Our formulation of the problem an d the numerical meth- 
ods ar e essentially the same as that of ITomimura & Eriguchil 
1 2005) , i.e. the num erically exact method, and that of 



IDuez & Mathij ( 120100 . i.e. the zero-flux-boundary method. 



2.1 Grad-Shafranov equation 

We calculate self-gravitating, axisymmetric, stationary mag- 
netized stars in order to obtain magnetized equilibria with 
strong toroidal fields in the Newtonian gravity. We assume that 
the system is in a stationary and axisymmetry state. For rotat- 
ing stars, the rotational axis and the magnetic axis coincide and 
the rotation is assumed to be rigid. The star has no meridional 
flows. The conductivity of the stellar matter is infinite, i.e. the 
ideal MHD approximation is employed. There is no magne- 
tosphere around the star. In other words, no electric current 
exists in the vacuum region. Therefore, the toroidal magnetic 
field is confined within the star and the only poloidal compo- 
nent can penetrate the surface and extend to the outside of the 
star. We use the polytropic equation of state as below: 

P = K P ^ N . (2) 

Here p, p, K and TV are the pressure, the mass density, the 
polytropic constant and the polytropic index, respectively. We 
fix TV = 1 for simplicity when we compute stationary config- 
urations by using the numerically exact method. This choice 
of TV is the same as t hat adopted in the previous works (e.g. 
lLander & Jo"nesll2009l) . 

We introduce the magnetic flux function "t in order to ob- 
tain magnetized equilibria efficiently. The magnetic flux func- 



tion is defined as below 
1 9* 



B 



B» 



i a* 



(3) 



r 2 sin 9 89 ' r sin 9 dr 

where B r and Be are the components of the magnetic field 
in the r-direction and the ^-direction, respectively. Using the 
flux function \l/, we derive the Grad-Shafranov equation from 
Maxwell equations as follows: 

% , d 2 t> sin0 d ( 1 <9*\ „ . „jV 



dr 2 



sin 9 



where j v and c are the ip-component of the current and the 
speed of light, respectively. 

The magnetic flux function ^ can be expressed as: 



* = r sin 9A a 



(5) 



where A v is the ^-component of the vector potential which 
satisfies V x A = B. The form of Grad-Shafranov equation 
can be rewritten as: 



47T 

A( A v sin ip) — j v 

c 



sm</p, 



(6) 



where, A denotes the ordinary Laplacian operator. This equa- 
tion seems strange because the system treated in this paper is 
axisymmetric and (,0-dependency would not appear. The rea- 
son why we introduce the term sin ip is that by introducing the 
sin 95 term we can 8 et an equation which contains the three 
dimensional Laplacian operator as shown above. The tree di- 
mensional Laplacian operator which emerges by this seem- 
ingly 'strange' device allows us to transform the equation into 
an integral form as shown below. By taking the boundary con- 
dition for the vector potential into account and using Green's 
function which satisfies the boundary condition, we derive the 
integral form of the GS equation as follows: 



ltp(r) smtp = 



j v (r') 



jr sin ip' d v + A V: h (r, 9) sin ip, (7) 



where A v .h is a homogeneous general solution to Eq.((6} as 
follows: 



A Vih {r,e) = Y, kt£rr + &»— ^(cos0). 



(8) 



Here, ro is a certain constant which is the stellar radius for 
spherical configurations, and a n and b n are constant coeffi- 
cients which are obtained by applying the boundary condi- 
tion. We will be able to obtain the stationary distributions of 
the magnetic vector potentials by solving this equation. Since 
these equations are of elliptic type partial differential equa- 
tions whose source term is j v , the boundary conditions are 
very important and have significant influences on the global 
structures of the vector potentials or the magnetic flux func- 
tions. We will deal with this problem about the boundary con- 
dition in Sec. [3] 



2.2 Toroidal magnetic fields 

Once we have obtained the flux function by solving the GS 
equation, it is easy to calculate the poloidal magnetic fields 
directly. On the other hand, we can obtain the toroidal com- 
ponent of the magnetic field by using the conserved quantity 
along the flux function which can be expressed by an arbitrary 
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function of "if. This arbitrary function appears because of the 
as sumption o f the axisymmetry. This function is called re(^fr) 
in lTomimura & Eriguchil (120051). F in lDuez & Mathij feOld) 
and T in IGlampedakis et al. (2012j ). In th is paper we name it 
k($) after Tomimura & Er iguchil ( 12005b . The toroidal com- 
ponent of the magnetic field is obtained from the following 
relation 



B« 



«(*) 



(9) 



This arbitrary function also appears in the expression of the 
current density as follows: 



J_dre(£) 

47T d* 



B + prsm9n( i S>)e if 



(10) 



where, //(\&) is another arbitrary f unction of ty. The arbi trary 
func tion n(^>) is the same as G in lDuez & Mathisl J2010f) and 
F in Glampedaki s et al.l J20121) . Then, we express the ip com- 
ponent of the current density as below: 



■Li 

c 



1 d«(*) «(*) 



+ prsm0(j,(ty) 



(11) 



47T d* r £ 

The first term of Eg. Jilt is the force-free current density part 
and the second term is the non-force-free current part because 
of V x B — 4:irj /c. If p = in a certain region, the magnetic 
field there is force-free, because j oc aB in that region. We 
will call the first term as the re term, and the second term as 
the p term of the current density in this paper. In a naive treat- 
ment, it seems to be enough to make the contribution from 
the re term larger in order to make the toroidal magnetic fields 
larger. However, it has been very difficult to make the influence 
from the re term strong. The reason for that is as follows. The 
distribution of ^ is obtained by solving the GS equation, but 
the source term of the GS equation contains the re term which 
is an arbitrary function of ^ and is confined to a restricted re- 
gion in the interior of the star because we impose the magnetic 
flux function to be smoothly connected on the stellar surface. 
Therefore, if we change the functional form and values of pa- 
rameters which appear in the functional form of re as well as 
the the region where re does not vanish, the distribution of the 
flux function also changes according to the changes of re. We 
will discuss this difficulty in Sec |4.2| 



3 SURFACE CURRENTS 

In this section, we deal with the relation between the surface 
current and the magnetic field, which is deeply influenced by 
the boundary conditions. The surface current can be defined 
either by the discontinuity of the derivative of the magnetic 
flux function, or by the homogeneous term in the integral rep- 
resentation for the vector potential. The both definitions for 
the surface current give exactly the same values as we will see 
later. 



3.1 Relation between the surface current and the 
discontinuity of the magnetic field 

At first, we will show a relation between the surface cur- 
rent and the discontinuity of the magnetic field as follows. If 







the Ampere's equation is applied to an area S bounded by a 
boundary dS, we can write it by an integral form as follows: 

<f B-d£=—[j-dS, (12) 

JdS c JS 

where dl and dS are a line element and a surface element, 
respectively. We apply this equation to an infinitely small area 
in the meridional plane of the star bounded by four lines as 
follows: 

Ar 

~T ' 

Ar 

e = 6>i, (is) 

e = 0i + A0, (16) 

where 9i, Ar and A9 are a constant, infinitesimal widths in 
the r-direction and in the ^-directions, respectively. For this 
infinitesimal area, we obtain 

[Bf x (r a ,e)-B i r n {r s ,e)JAr+ \Bf{r a ,6) - Bf-(r s ,8)\r 3 A8 

47T /•'-a+Ar/2 



(13) 
(14) 



Jr,-A? 



/2 



j v (r,6)rdrA0, 



where B ex and B vn are the exterior and interior values 
of the magnetic fields. Here, if the current density is a sur- 
face current on the stellar surface (r = r s {9)) defined by 
j v (r, 9) = j S ur(0)8(r — r s (6)), we can integrate the equa- 
tion as below 



r,+Ar/2 

s -Ar/2 
j's+Ar/2 

s -Ar/2 



j v (r,9)rdr 



jsur(0)5{r- r 3 )rdr = r s j aur (9). (18) 



Since the r dependence of the magnetic fields are continuous 
on the stellar surface (r = r 3 (9)), we obtain a relation as fol- 
lows: 



B("(r.,J)-B, m (r.,l 



47T 

J sur I v I . 
C 



(19) 



The surface current is expressed by the discontinuity of the 9- 
component of the magnetic field. For more general situations, 
we can obtain the following equation for the surface current 
using the parallel component B\ | to the stellar surface, 



47T 



Bff{r a ,0) - B*"(r.,0) = — j mr {6) 



(20) 



If the surface current exits, the parallel component of the mag- 
netic field must be discontinuous. We emphasize that the value 
of the discontinuity of the magnetic field between just inside 
and just outside of the stellar surface equals the surface current 

density. 

iGlampedakis et al.l J2012T) expressed the surface current 
density in a different way as follows. They defined the surface 
current by imposing the discontinuity of poloidal magnetic 
fields at the stellar surface (see Eq.67 in IGlampedakis et al.l 
|2012|) . Their discontinuous boundary condition is just an as- 
sumption without a firm foundation as follows: 

bf = m X ^ bT - bf = (1 - iW = jsur. (21) 
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Figure 1. Poloidal magnetic field lines with a purely dipole (n=l) surface current (left) and a purely quadrupole (n = 2) surface current (right). 



Here b and £ indicate the magnetic field and a discontinu- 
ity parameter, respectively, in their paper. Since their model 
is a purely dipole configuration, the exterior solu tion of bg i s 
b e e x (x) = ± da/dx sinfl (Eq.61 in lGlampedakis et alj|2012h . 
Here a and x indicate the r-component of the flux function 
and the dimensionless radius normalized by the stellar radius. 
Then we can calculate the distribution of their surface current 
density as follows, 



(1 



. 1 dot . 
^dx-^ 



-jo sin 9, 



(22) 



(s=i) 



where 



Jo = 



(£-1) 



J_d& 
2x dx 



(x = l) 



Since they calculated only models with £ > 1 in their pa- 
per, the surface current density of their models flows in the 
opposite direction to the interior bulk toroidal current density 
inside the star. In other words their surface current is nega- 
tively flow ing with respect to the bu lk interior currents. Ac- 
cording to lGlampedakis et alj I I2012I) . if the value for the dis- 
continuity for the poloidal magnetic field is increased, the en- 
ergy of the toroidal magnetic fields becomes larger (see Fig. 
5 in Glampedakis et al. 2012). Therefore we conclude that the 
negative surface current sustain strong toroidal magnetic fields 
comparing with those in Tomimura & Eriguchi models with- 
out surface currents. 



We can obtain the magnetic field by calculating Eq. 0, 

A ■ 1 fjsur(r',0') . , ,z l 

A v sin (p — - I — ; — sin (par. 

c. \r-r'\ 



(23) 



Since the surface current exists on the stellar surface r = r 3 , 
we can describe the surface current density using the Dirac's 
delta function: 



jsur(r',e') 



5{r' - r s )j sur {6') , 



(24) 



where j 3ur (9) is the surface current which flows along the sur- 
face. We can expand and integrate Eq,ll23t using the Legendre 
functions and the axisymmetry of the system. We obtain the 



AAr,0) = J2 



2tt 



-Pn(cos6)f„(r,r s )r 2 3 



—J n(n + 1) 

sin e'P^{cos 9') j sur (9')d9', (25) 



(26) 



where /„ is a function defined by 

f(rr)-t r?/r " +1 ' (r " rs) 
fn(r,r s )- < r „ /r n+i^ ( r < rs ) 



Now we calculate the vector potential and the magnetic field 
by giving a ^-distribution for the surface current density. We 
assume the ^-distribution of the surface current can be ex- 
pressed by the expansion using Legendre functions, 



3.2 Surface currents in the integral representation 

Using the integral representation for the vector potential, we 
can see the surface current from a different point of view. We 
assume that a magnetized star has no currents in the stellar 
interior except for the surface current. It implies that the source 
term for the GS equation consists only of the surface current. 



, . ^— \ In +1 „i , 

3sur{9) = 2_^ — a„P n (cost 

n = l 



(27) 



where Q n 's are dimensionless coefficients related to the nth 
associate Legendre function P^(cosf?). Then using the or- 
thogonality among the Legendre functions, 



,(n + l)n 



/ sin9' PUcos9')Pl(cos9') d9' = 2- 
Jo 2n + 1 



(28) 
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we can obtain an analytically expressed solutions as follows: 



and 



A v (r,0) = y^ ot n rlf n (r, r s )P*(cos I 



and 



*(r, 6) — ^ a„rr 2 s f n (r, r s )P*(cos9) sin6>. 



(29) 



(30) 



If we set a n = a n + 6 n , the right hand side of Eq. J29t 
is exactly the same as homogeneous general solutions of Eq. 
{8]l. Therefore, adding the homogeneous term to the inhomo- 
geneous solution of the GS equation corresponds exactly to 
adding the surface current at the boundary surface. 

Fig. [U shows the poloidal magnetic field lines for config- 
urations with the purely dipole (n = 1) surface current (left 
panel) and with the purely quadrupole (n = 2) surface current 
(right panel). It should be noted that each model has no interior 
currents except for the surface currents as we have described 
in this section. 



4 NUMERICALLY EXACT CONFIGURATIONS FOR 
OPEN MAGNETIC FIELDS WITH SURFACE 
CURRENTS 

We will show numerically exact structures of magnetized stars 
with open field lines and with/without surface currents in this 
section. At first, we will display configurations which have no 
surface currents. Although they are the same as those obtained 
in lYoshida & Erig uchi (2006) and [Lander & Jones (2009), we 
will check these models from a different point of view, i.e. 
in the context of the influence of the surface currents. Next, 
magnetized stars which have surface currents will be treated 
and discussed. 



4.1 Setting of the problem 

As discussed before, we have solved the integral equation de- 
rived from the GS equation by considering the presence of the 
surface currents as follows: 



A v (r) simp — — I 



1 f j v (r') + jsur(r') , , , 



\r — r'\ 



sin(p d r , (31) 



where j sur is the surface current density of the magnetized 
star. We choose the following two different distributions for 
the surface currents: 



jsur{r,0) 

c 
and 
j S ur(r,0) 



-josm85(r — r s (0)), (dipole distribution) , (32) 



As for the arbitrary functions appearing in our formula- 
tion, we choose the following forms: 



/ 



M(*) = Mo , 



(34) 
(35) 



«(*) 



0, 

KQ 
fc + 1 



for* < *Vmax, 
(f-^Vmax)^ 1 , for* > * Vmax , 



K ' W ^ K0 (*-*V maX ) fc 



for* < * ymax 
for * > * 



V max j 



(37) 



Here, fio, k and kq are constant parameters and *v max means 
the maximum value of * in the vacuum region. In this paper, 
we fix k — 0.1. These functional forms a nd the value of k are 
the same as those chosen in other papers (lYoshida & Eriguchil 
|2006J; lLander & Jonesll2009l). In this section, w e set the poly- 
tropic index N — 1 (e.g. lLander & Jonesll2009D . and q = 0.9 
where q is the ratio of the polar radius to the equatorial radius 
defined by (see F ujisawa et alj2012l) 



r s (6 = 0) 
r 3 ((9 = tt/2)' 



(38) 



-JO sin9cos85(r — r s (8)). (quadrupole distribution) (33) 



Concerning the angular velocity fi, we choose values of fi = 
constant = Clo — 1.0E-2 for rigidly rotating configurations 
and Q,q — 0.0 for non-rotating models. Here quantities with" 
represent the corresponding ones transformed i nto dimension- 
less f orms as shown in Appendix IA2I (see also lFuiisawa et al.l 
120121) . 

Although the equations of st ate influence the strengt hs of 
the toroidal magnetic fields (see lKiuchi & Kotakdl2008i ). we 
have treated only N — 1 polytropes because the main concern 
in this paper is how the surface current density affects the dis- 
tributions of the magnetic fields, in particular to the toroidal 
magnetic fields. In order to examine the accuracies of solu- 
tions, we have used the virial relation as shown in Appendix 
[A3] 



4.2 Configurations without surface currents 

Since the value of ko affects the local behavior of the toroidal 
magnetic field di stributions, in particular on its maximum 
value (see lLander & Jo nes 2009), we have solved the magne- 
tized configurations by changing the value of ko for two values 
of fio- 

Obtained results are plotted in Figf2]which shows the ra- 
tio of Mp/M against ko. As seen from this figure, there is a 
minimum value of the ratio. It implies that the toroidal mag- 
netic field energy increases to its maximum value at ko ~ 40 
for non-rotating configurations and at ko ~ 45 for rigidly ro- 
tating models. Since the term related to the rotation does not 
appear in the current density formula (see EqllOl. the rotation 
affects the toroidal magnetic field distributions only slightly. 
Therefore, we will display and discuss only configurations 
with rotation in the following part of this section. 

In many investigations which have been done by apply- 
ing numerically exact methods or by structure separated GS 
solving method zero-flux-boundary methods, almost similar 



result s as ours shown in Fi g. |2]have b een obtained (see Table 
2 in lLander & Jonesll2009l F ig. 12 in lCiolfi et af1l2009l. Fig. 
4 in Glampedakis et al. 2012). Therefore, this behavior of the 
toroidal magnetic field against the value of /to is likely to be 
a general feature of stationary magnetized stars which have 
open magnetic fields. 

In order to consider the reason of the presence of these 
minimum values, in Fig[3] shown are the distributions of the 
magnetic field components (left panels) and those of the two 
(361 com P onen t s of the current density formula (right panels). Dif- 
ferent curves in the left panels mean B r (9 — 0) (dotted line), 
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Figure 2. The ratio of the poloidal magnetic energy to the total magnetic energy against the parameter «o ■ The solid line denotes a rigidly rotating 
sequence with Qq = 1.0E-2.0 and the dashed line denotes a non-rotating sequence, i.e. a sequence with f2 = 0.0. Minimum values seem to be 
attained at kq ~ 40 (for the non-rotating sequence) and kq ~ 45 (for the rigidly rotating sequence). 
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Figure 3. Left panels: Distributions of the components of the magnetic field [B r (8 = 0), B$(0 = 7r/2), B v (9 = ?r/2)] are shown. Right panels 
: Distributions of the two source terms on the equator, which appear in the current density formula, are shown. Chosen values of k are ko = 10 
(top) and kq = 65 (bottom). 
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B 9 (9 = 7r/2)(dashed line) and B^{6 = 7r/2)(solid line) dis- 
tributions for the region f [0 : 1.0]. In the right panels, the solid 
line denotes the force-free k term and the dashed line denotes 
the non-force-free fi term in the current density formula. 

As seen from these panels, by increasing the value of kq, 
which corresponds to increasing the maximum strength of the 
toroidal magnetic field, from top panels to bottom panels, we 
can find that the width of the toroidal magnetic field region be- 
comes narrower. The values of the toroidal magnetic field en- 
ergy seem to depend on the distributions of the toroidal mag- 
netic fields and of the current densities, in particular, on the 
maximum value and the width of the toroidal magnetic field 
distribution. Although the maximum value of the n term of the 
current density increases with ko, its width becomes narrower 
as ko is increasing. It implies that the slope of the distribution 
of the k term becomes steeper for the large value of the cur- 
rent density. This tendency is also seen in the distribution of 
B v . By contrast with this, the distributions of the /j, term are 
almost the same because we fix q and Clo which are related to 
the characteristic nature of the non-force-free magnetic fields. 

From Eq. dlOt and our numerical results, we can find that 
in order to sustain strong toroidal magnetic fields (appearing 
in the right hand side of the equation), the strong toroidal cur- 
rent density (appearing in the left hand side of the equation) 
is required. It should be noted that the strength of the toroidal 
magnetic field seems to be related deeply to the strength of the 
total current density. This can be seen from the following ar- 
gument. We introduce several definitions about integrated cur- 
rents as follows: 



(+) 



■h 



f(-) 

J tot 



f K 

•Jtot 



J? - 

°tnf 



J ( v +) dS , 



j ( v ~ ] dS. 



j*dS . 



]ys , 



(39) 
(40) 
(41) 
(42) 



Jsur = / jsurdS — -jo r s (9)sin9d9, (43) 
J s mcr Jo 

Jtot = Jtot "t" Jtot + Jsur — Jtot + Jtot + Jsur ,(44) 

where S mer denotes the meridional plane which is perpen- 
dicular to the (^-coordinate and dS is an area element in 
the meridional plane. Here, j^ , j^~ , j'S and J^ the tp- 
component of the positively flowing interior current density, 
the (^-component of the negatively flowing interior current 
density, the n term of the current density and the /i term. Fur- 
thermore, j\rf, Jtot , Jtot, Jtof Jsur and Jtot are the to- 
tal positive bulk interior current, the total negative bulk inte- 
rior current, the total n term bulk interior current, the total /i 
term bulk interior current, the total surface current and the total 
(bulk + surface) current in the meridional plane, respectively. 
As we shall see, these quantities will play key roles to under- 
stand the problem. 

In Fig[4] the total current, the total fi current and the to- 
tal k current of the star is plotted against ko- We find from 
Fig. E] that the total current does not increase as «o increases 
and the total k current increases to its maximum value near at 
ko ~ 45. We will denote the maximum value of the total cur- 



rent as J t ^ in this paper. This ko value is the same as that 
for the minimum value of M v j ' M. Therefore it is important 
to note that there is an upper bound of the total current for con- 
figurations if we consider a stationary sequence with different 
values of ko • 

This upper bound comes from our boundary condition for 
the current density. Since we have imposed that the outside of 
the star is vacuum, the current density needs to vanish in that 
vacuum region. As we have seen, the magnetized stars need 
large and strong toroidal currents in order to sustain strong 
toroidal magnetic fields. However, the boundary condition sets 
limit to the total current of the star as seen from Fig|4] As 
a result, the region where the current density attains a rather 
large value becomes smaller and the slope of the distribution 
of the current density becomes steeper in order to sustain the 
stronger toroidal magnetic field in the narrower region. 

Moreover, larger values of ko cause the maximum value 
of the magnetic flux function in the vacuum region larger, in 
general. As far as our boundary condition for the magnetic 
flux function to be smooth at the stellar surface is employed, 
the support of the k function becomes smaller and smaller as 
the value of ko is increased. In other words, increasing the 
value of ko might, in ordinary situations, result in increasing 
the interior currents but at the same time decreasing the sup- 
port region of the function k because the maximum value of 
the magnetic flux function in the vacuum region also becomes 
larger as explained before. 

This is the reason why in the present investigation as well 
as in other works thus far done nobody could obtain solutions 
which exceed this upper bound. To overcome this limitation 
about the size of the confined region of the large toroidal mag- 
netic field, the magnetized stars needs other kinds of distribu- 
tions for the toroidal current densities. 

From these consideration, we need to devise some means 
to fulfill the following seemingly contradicted requirements at 
the same time. 

(1) (^-currents must be increased. In ordinary situations, 
this would results in reducing the support region of the func- 
tion k because of the increase of the maximum value of the 
magnetic flux function in the vacuum. (2) The support region 
for the function must be widened. In ordinary situations, the 
support region of the function k is wider for the smaller values 
of Kq. 

These two seemingly contradictory requirements could 
be realized by introducing negatively flowing currents near/on 
the surface because the negatively flowing currents allow the 
positively flowing interior currents to become larger and at 
the same time negatively flowing currents near/on the surface 
could reduce the value of the magnetic flux function in the 
vacuum region and result in the smaller value for ^v-max- 



4.3 Configurations with surface currents - Dipole 
currents 

As explained in the previous subsection, in order to exceed 
the upper bound of Jtot found in this paper and to reduce the 
value of tyvmax, we will try to investigate the magnetized 
stars which contain oppositely flowing surface toroidal cur- 
rents against the interior 'bulk' currents which are flowing in 
a certain direction. We will call such oppositely flowing cur- 
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Figure 4. Left panels: The total current Jtot (thick solid line), the total fi term interior current Jt ot (dashed line) and the total K term interior 
current J t K ot (thin dashed line) are plotted against kq for configurations without surface currents. Right panels: Only the J? L is plotted against kq. 
The magnitude of the J^ ot has an upper bound near at kq ~ 45 . 
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Figure 5. The ratio of the poloidal magnetic energy to the total magnetic energy of the models with surface current against the parameter jo . 
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© 0000 RAS, MNRAS 000, 000-000 



10 K. Fujisawa & Y. Eriguchi 



0.2 

0.15 

0.1 

0.05 



-0.05 

-0.1 

-0.15 

-0.2 



jo = 2.5E-3 



-"phi 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

r 

jo = 7.5E-3 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

r 




0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




jo = 7.5E-3 





■ , ' 


/ \ K 

























0.2 



0.4 



0.6 



0.8 



Figure 7. The same as Fig:3 except for configurations with non vanishing surface currents, i.e. configurations with jo = 2.5 X 10 3 (top) and 
jo = 7.5 X 10~ 3 (bottom). 



rents as negative currents, hereafter. In short, we will assume 
that there could be toroidal surface currents which flow to the 
negative direction compared to the flow direction of the inte- 
rior main currents which we will call the interior bulk currents. 
In fact the effects of the presence of the oppositely flowing 
surface currents are similar to those of the boundary condi- 
tions treated in section 5.2 of iGlampedakis et alj 1 120 121) as we 
have shown in Sec l3. II As we shall show, magnetized stars 
with negative surface currents will be able to sustain much 
stronger interior bulk currents and have much stronger toroidal 
magnetic fields because the oppositely flowing surface current 
cancels the effects of the interior toroidal currents to certain 
extent and results in configurations which have the following 
two special characteristics. 



(i) In such configurations, although the total currents Jtot 
are below the upper bound J^ ax discussed before, much 
stronger positive interior bulk currents J tot are allowed to ex- 
ist. 

(ii) At the same time, in such configurations, the absolute 
values of the magnetic flux function in the outer vacuum re- 
gion can become smaller than those of configurations without 
negatively flowing surface currents. Thus the support region 
for the arbitrary function k^) can become wider than that for 
configurations without surface currents. 



In this subsection, as an example, the dipole-like distri- 
bution for the surface current as Eg. (132) is employed. If the 
magnetized stars are purely spherical with no interior cur- 
rents within the stars, dipole-like surface currents result in uni- 
formly distributed interior magnetic fields and purely dipole 
exterior magnetic fields (see FigQ}. Thus, if the surface current 
densities are much stronger than the interior current densities, 
the interior magnetic fields become almost uniform and there 
are no closed magnetic fields inside of the stars. The toroidal 
magnetic fields could not appear in such configurations. 

We have calculated many stationary configurations with 
surface currents for several values of ko . As a result, we found 
that a model with ko = 10 has the smallest value of M p /M 
in all our stationary solutions. Thus we will show only con- 
figurations with ko = 10 in this paper, but it should be noted 
that all other models show almost the same tendency as that 
of configurations with ftn = 10 which we will describe below. 
Fig.|5]shows the values of M p /M and Fig.|6]shows the val- 
ues of Jtot (thick solid line), J tot (thin solid line) and J sur 
(thick dashed line) in the left panel and the value of J tot (thick 
solid line), J^ ot (thin dashed line) and Jf ot (thin solid line) in 
the right panel against the values of jo for configurations with 
ko = 10. In these models, there is no negative current j^' in 



the star. Therefore, j[ ot 



and J, 



(+) 



Jtot + Jtot in theSe 



configurations. From this figures, we can see that if we in- 
crease the value of jo, the total bulk current of the magnetized 
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star, J tot , becomes larger (Fig(6) and the value of M p /M 
becomes smaller (Fig(5}- We can see from left panel of Fig|6] 
that the total current J tot is almost the same as the upper bound 
of the total current J t ™ defined before. However, the to- 
tal positive current J toi becomes much larger than this upper 
bound. Especially noted from the right panel in Fig(6] the total 
K current term becomes much larger and the total /j, current 
term becomes slightly small. This can be considered as an ev- 
idence that the negative surface current cancels some contri- 
butions of the positively flowing interior toroidal current from 
the k term current. It is remarkable that the value of M p /M 
attains about 0.7 when jo — 7.5 x 1CP 3 . It implies that the 
magnitude of the toroidal magnetic field energy is almost the 
same order as that of the poloidal magnetic energy for those 
models around j — 7.5 x 10" 3 . The ratio M p /M ~ 0.65 
seems to be the minimum value in the present parameter space 
because we could not succeed in obtaining numerical solu- 
tions when jo > 1-0 x 10~ 2 . Since the surface current with 
jo = 2.0 x 10~ 2 should be considered tremendously strong, 
their fields would become nearly uniform in the interior and 
purely dipole in the outside of the star by the surface current 
as seen FigQ] Moreover, when there are no closed poloidal 
magnetic field lines inside the stars, the magnetized stars can- 
not have toroidal magnetic fields because the poloidal mag- 
netic fields are originated from the closed current densities 
which are assumed to be parallel to the closed poloidal mag- 
netic fields as seen from the current density formula. 

Fig.|7]shows the components of the magnetic fields (left) 
and the non-force-free and the force-free term in the cur- 
rent density formula (right) are plotted against f. We choose 
jo = 2.5 x 10~ 3 (top panels) and ] - 7.5 x 10" 3 (bot- 
tom panels). From these panels, we can see both the width and 
the strength of the k term are increasing as jo, but that the 
strength of the /i term does not change so much. This result 
means that the oppositely flowing surface current affects only 
the force-free term in the current density formula significantly. 
As discussed before, since the /j, term is non-force-free and af- 
fected mainly from the global structures of the stars, i.e. by the 
value of the axis ratio q. The distributions of /i term in the cur- 
rent density formula would not show large change. We have 
computed these configurations by fixing the values of q and 
Qo and the total magnetic strengths of the stars are nearly the 
same for the different configurations. Therefore, the distribu- 
tions of the fi terms are nearly the same for different values of 
Jo- 
lt is needless to say that in order to increase the total cur- 
rent keeping the /i term nearly the same, the k term must be- 
come larger and stronger as we see in these panels and Fig|6] 
Since the stronger and steeper distributions of the k term re- 
sult in the stronger toroidal magnetic fields, the presence of the 
oppositely flowing surface current should be required for the 
larger and stronger toroidal magnetic fields. In other words, 
the oppositely flowing surface current density can sustain the 
strong toroidal magnetic fields. The maximum value of the 
toroidal magnetic field and the size of the region where the 
most of the toroidal magnetic field exists can be increased by 
adding and increasing the surface current densities. 

In Fig[8]the distributions of the magnetic flux functions 
on the equator (8 = 7r/2) for two configurations, one without 
surface current (left panel) and the other with surface current 



(right panel), are displayed. As seen from left panel, the value 
of the magnetic flux function at the equatorial surface for the 
configuration without surface current becomes bigger as the 
ko increases, because the maximum value of the flux function 
becomes bigger. On the other hand, the value of the magnetic 
flux at the equatorial plane with surface current (right panel) 
does not change very much even if the values of *l/ m ax in- 
crease. Therefore, the negative surface current make the flux 
function at the equatorial surface smaller than that for the con- 
figuration without surface current. This reduction of the value 
of the magnetic flux function in the vacuum region can allow 
the wider support region to exist. 

As for the geometry of magnetic fields, the surface cur- 
rents bend the poloidal magnetic fields on the stellar surface 
as we have described in Sec l3.ll Fig(9]shows that the poloidal 
magnetic field lines bend due to the presence of the surface 
currents for configurations with jo = 2.5 x 10 -3 (left) and 
jo = 7.5 x 10" 3 (right). 

Thus far, nobody has obtained configurati ons with the 
surfac e currents in previous works (see Fig. 2 in Fujisawa et al. 
120121) . On the other hand, configurations shown in Fig|9]have 
discontinuities on their surfaces. The directions of the discon- 
tinuities depend on the directions of the surface currents. In 
the northern hemisphere, the outward poloidal field lines are 
bended to the left side by the oppositely flowing surface cur- 
rents. On the other hand, if we add the surface currents whose 
flowing direction is the same as that of the interior bulk cur- 
rent, they are bended to the right side and the toroidal mag- 
netic fields become weak because the surface currents flowing 
to the same direction as the interior bulk currents work so as 
to reduce the strengths of the interior bulk currents. The dis- 
continuitie s of the magnetic fields in our models are the same 
as those of lGlampedakis et alj ( 120120 (see Fig.5 in their paper 
carefully). In fact, what they did in their paper, i.e. by impos- 
ing structures in which the magnetic fields on the stellar sur- 
faces have some amount of discontinuities are effectively the 
same as adding the oppositely flowing surface currents to the 
magnetized stars. 



4.4 Configurations with surface currents - Quadrupole 
currents 

We consider configurations with other surface current distri- 
butions. We add the surface currents expressed by Eq.(|33t, 
which results in the quadrupole distribution of the poloidal 
magnetic fields. However, the toroidal magnetic field for this 
surface current cannot become large enough. As the strength 
of the surface current is increased, we get configurations with 
no toroidal magnetic fields whose poloidal magnetic fields are 
not closed inside of the stars (see the right panel in Fig[T}. The 
toroidal component of the magnetic field vanishes in such a 
configuration. 

In order to sustain strong toroidal magnetic fields, we 
need strong toroidal currents in the stellar interior as discussed 
for the dipole-like distributions of the surface currents. How- 
ever, as seen from the results for the quadrupole-like surface 
currents, surface currents contain both the negative compo- 
nent and the positive component in the surface currents. More- 
over, the strengths of those oppositely flowing currents are the 
same. Therefore, the total surface current due to the purely 
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Figure 8. The distributions of magnetic flux function on the equator (9 = 7r/2). Left panel: These lines denote k = 0.0 (solid line), ho = 20.0 
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Figure 9. Distributions of the poloidal magnetic fields. Parallel components of the poloidal magnetic fields along the surfaces are discontinuous on 
the stellar surfaces. The values of M P /M ~ 0.8785 (left) and ~ 0.6979 (right). 



quadrupole surface current of the purely spherical star van- 
ishes as 



L 



Jsur@>d 



r s jo sm a cos t 



= 0. 



(45) 



Consequently, this kind of surface current cannot counteract or 
cancel the effect of the interior bulk current density. It implies 
that we need 2 e -th poles with I = odd numbers of the magnetic 
moments (dipole, octopole etc.) or locally strong surface cur- 
rents which are not anti-symmetric about the equatorial plane 
in order to sustain the strong toroidal fields. 



5 REASONS FOR APPEARANCE OF STRONG 
TOROIDAL MAGNETIC FIELDS — 
COEXISTENCE OF OPPOSITELY FLOWING 
MULTI-v>CURRENTS 

In the previous section we have discussed that the presence 
of negative (oppositely flowing) surface currents in addition 
to positive interior bulk currents could allow more interior 
currents to exist within the stars. In particular, large interior 



toroidal currents could be realized by introducing negatively 
flowing surface currents in addition to positively flowing inte- 
rior currents. Consequently, such negatively flowing toroidal 
currents lead to larger positively flowing total currents, J^„ t , 
within the stars, although the total currents, Jtot, have their 
upper bound as explained before. Thus it is these larger pos- 
itively flowing interior currents, J t( J , that cause the toroidal 
magnetic fields stronger. 

Among stationary magnetized stars thus far obtained in 
many papers, some authors have found configurations with 
large toroidal magnetic fields by treating the problem dif- 
ferently (seefloka & Sasaki 2004, Duez & Mathis 20j3 and 
lYoshida et alj2012l) . However, no authors have explained rea- 
sons why there can appear such magnetic fields with large val- 
ues of the toroidal magnetic field energies. 



5.1 Zero-flux-boundary approach: N = magnetized 
spherical configurations 

In this subsection, we will reconsider simple configurations 
with large toroidal magnetic fields from our standpoint of tak- 
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ing the negatively flowing currents into account. In order to 
clarify the reasons for existence of large toroidal magnetic 
fields, it would be helpful to employ as simple configurations 
as possible. 

As for the mechanical structures of the magnetized 
stars we consider TV = polytropes, i.e, incompress- 
ible fluids. Although the magnetic fields might become very 
strong, the shapes of the stellar configurations are assumed 
to be spheres. It is possible to consider that strongly mag- 
netized stars have spherical surfaces by confining all the 
magnetic fields within the surfaces. This situation could 
be realized by treating the closed magnetic fields not only 
for the toroidal fields but also for the poloidal fields. Un- 
der these assumptions, the formulations used by several au- 
thors ( Chandrasekhar & Prendergast 1956; Prendergast 1956; 
Duez & Mathis 2010; Glampedakis et al. 2012; Yoshida et al. 
20121) can be applied. Since our purpose of this part of the pa- 
per is not to obtain 'new' stationary configurations but to un- 
derstand the reasons for appearance of large toroidal magnetic 
fie lds, we will follow mostly th e zero-flux-boundary sc heme 
of lPrendergastl i 19561) (see also iDuez & Mathisl J2010l) ). but 
we take account of the presence of negatively flowing surface 
currents as well as negatively flowing interior currents in ad- 
dition to positively flowing interior currents. 

As for the arbitrary functions, we choose the functional 
form of ft as follows: 



ft(*) = K *, 

and consequently 



dft(*) 



fto- 



(46) 



(47) 



This choice for ft^) as well as the form for /i(^E') = 
flo (constant) have been commonly used in almost all zero- 
flux-boundary approaches which have tre ated configurations 
with closed dipole magnetic fields (see 



Ioka & Sasaki 2004, Duez & Mathis 2010 



Prendergast 1956, 
and I Yoshida et alj 



2012). From Eq.dlOt. the form of the current density becomes 
as below: 



c 



■2 

'mi 



4-7rr sin 9 



ty + pr sin 9 fio . 



(48) 



Since the GS equation is an elliptic type partial differ- 
ential equation of the second order, we need to impose bound- 
ary conditions to obtain solutions consistently. We should note 
that, in all zero-flux-boundary approaches thus far carried out, 
the constant fto plays a role as an eigenvalue of the problem be- 
cause boundary conditions have been imposed at finite places 
in the space in most investigations. One example of the bound- 
ary conditions might be as follows: 



*!r=0 = 



*l 



, (at the centre), 

, (on the stellar surface) 



(49) 
(50) 



» 



where r s is the radius of the spherical incompressible mag- 
netized stars treated in this section. It should be noted that so- 
lutions with d\&/dr| _ («) = in our models, which can be 
found only after we have obtained stationary configurations 
and checked values of the derivative d"$> /dr\ _ ( s ) for all the 
models, are essentially the same as those of IDuez & Mathisl 
(2010). It should be stressed once again that solutions which 



satisfy the condition d& /dr\ _ ( s > =0 would not always be 
found. It would be fortunate if one could find such solutions 
not by imposing that condition as one of boundary conditions 
but by just calculating solutions with the boundary conditions 
g9> and <50l 



5.2 Magnetic flux functions for spherical incompressible 
fluids with magnetic fields confined within the stellar 
surfaces 

In this section, we continue t o follow mostly the fo rmulation of 
IPrendergasJ dl956l) (see also lDuez & Mathisll2oT(il) . in which 
the surface currents were not taken into account explicitly, but 
in this paper we include the surface currents explicitly by mod- 
ifying their formulation. 

They treated incompressible fluid stars, i.e. N — poly- 
tropes by specifying arbitrary functions as we have already ex- 
plained before. Although incompressible stars seem far from 
realistic situations, from the standpoint of considering oppo- 
sitely flowing currents including surface currents stressed in 
this paper, it is very useful to be able to get such solutions and 
discuss the role of the oppositely flowing currents analytically. 
Nevertheless, in this paper, we will also compute TV = 1 poly- 
tropes numerically and discuss the effect of the equation of 
state. 

For the functional forms we have chosen, the Grad- 
Shafranov equation, Eq.©, can be written as below: 

A** + fto* = -4ivfi pr 2 sin 2 9 = S(r, 9). (51) 

where p is the averaged value of the density. It should be noted 
that this is a linear equation for the magnetic flux function *. 
When p is constant throughout the stellar interior, we can inte- 
grate this GS equation easily by expressing the solution in the 
integral form by using Green's function and spherical Bessel 
functions and Gegenbauer polynomials as follows: 



1=0 

= — 47T^0 2j K ^,l 
1=0 



(52) 



2^ + 3 



2(1 +!)(£ + 2) 



sin 9C e (cos ( 



W*+l I K ,t-j^ 

' s \ ' s 



r'Y e+1 [ ko,£— J pdr' 



+ — y W ( K( V— ) / r ' J W ( *o.«-w J Pdr' 



S(r' ,cos9')C e (cos9')dcos9' , 



*» = E 



r / r 

Ki/Ko,e—nsJe+i [ K o,t—rpr 
r« \ re 



T I T 

+ K2,eKo ] i—-rYt + i I Ko,e—nx 
ri ' \ r 3 



sin 9C t (cos ( 



Here * s and ^ denote the inhomogeneous solution and the 
homogeneous solution to the GS equation, respectively and 
Ji 's and Ye's are the spherical Bessel functions of the first kind 
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and the second kind, respectively and C e 's are the Gegen- 
bauer polynomials, ^i, Kq,i and K\ t i, Ki,t denote the £-th 
expansion terms of the magnetic flux function, the eigenvalues 
corresponding to k for the £-th component equations appear- 
ing in the expansion of the magnetic flux function as above, 
and coefficients of homogeneous solutions, respectively. It 
should be noted that, here in this paper, we consider only the 
dipole term (£ — 0) which can be considered as the sim- 
plest configuration for the spherical incompressible magne- 
tized fluid star. Moreover it is important to note that even such 
simple configurations contain the essential natures of the con- 
figurations we are seeking to understand. 

As we have described before, we impose boundary con- 
ditions for vp. One condition is >I> = at the centre of the star. 
This condition can be fulfilled simply in our situation here by 
setting K% — 0, because the spherical Bessel function of the 
second kind Y( = o does not vanish at the centre (r = 0) for 
the homogeneous solutions. As a result, we obtain the general 
expression of the solution as follows: 



will write its constant value as jo 
1 d#(r,0) 



4-7rr sin 9 dr 



Jo- 



(58) 



Thus the distribution of the surface current can be written as 
below: 



Jsur 
C 



- Jo sin f 



(59) 



* = K^Ko—r^Jl ( Kq 



» 



sin 9 



47r/ioKo sin 9 



,(•> 



Jl KQ 



r 



Vi Uo 



r 



- 13 

pr 



dr 



rl 



r 



K 0-r-T / Jl k 0-t-t 



r 



- '3 j ' 

pr dr 



Explicit forms of Ji, Yi and C are as below: 



3/2 



Ji( A ) = T^( sln -^" 
A^ 



A cos A), 



Yi(A) = — — ^-(cos A + A sin A), 
A 



-,3/2 



Cq' (cos 6) = 1. 



(53) 



(54) 



(55) 



We denote Kq as K o and K° as K\ for simplicity in the fol- 
lowing part of this paper. 

Next, we impose the other boundary condition '3/ = at 
the stellar surface. This condition is written as follows 



K\ Ji (kq) = 4-irfioYi (k ) 



Jl Kq 



» 



pr ni dr'. (56) 



From this equation we can obtain a relation between ko and 
K\ of our problem at hand. Thus just by giving either ko or 
A'i, one complete solution in our problem can be obtained. 
This is a nice feature of the simplest configurations, i.e. the 
N — polytropic configurations only with the £ = compo- 
nent for the confined poloidal closed magnetic fields. 

Finally, we will derive the surface current for our prob- 
lem. The homogeneous term of this solution is related to the 
surface current as we have calculated in Sec 13. 21 Thus the sur- 
face cmi'ent is expressed as 



Explicit forms of K\ and jo can be found in Appendix 151 

Therefore, if the star has negative surface currents, values 
of the magnetic flux functions in the large part of the stellar 
interiors are positive because of d^ /dr < on the stellar 
surfaces. 

As explained before, we can calculate one eigenfunction 
just by choosing one value of ko. By changing the value of 
ko and calculating the corresponding solution for \&, we have 
a series of solutions which are shown in Fig[l0]and Figllll 
FigQUdisplays the various total currents Jtot- 

Fig. Qo] shows how jo behaves for different value of Ko 
(left) and how A4 P /A4 behaves against the value of ko (right). 
We set the parameters satisfying the relation r^ 3 ' flop — — 1, 
,As seen from these panels, the function j'o(ko) has two spe- 
cial solutions which contain no surface currents at ko ~ 5.76, 
and at k ~ 9.10. Furthermore, there appear two singulari- 
ties along this curve at ko ~ 4.49, and at kq ~ 7.73. The 
sign of the magnetic flux function changes at the singularities. 
The value of M p /M decreases as the value of ko increases 
until it reaches the first solution without a surface current at 
ko ~ 5.76 and the value of the A4 P /A4 reaches its minimum 
of M p /M — 0.417. Hereafter we will call the eigenvalue kq 
of the first solution without a surface current as kJ„ and the 
eigenvalue ko of the second solution without a surface current 
as K^. We also denote values of Ko at the first and the second 
singular points along this curve as K* and k 2 s , respectively. 

We will discuss the behaviors of \l/ and jo only around 
K X m in this paper. However we see almost the similar be- 
havior s for kq ~ K 2 m , which one can find in several pa- 
pers jlloka & Sasakil2004llDuez & Mathisl2010l,lYoshida et al.1 
l2012h . 

We can classify the eigen solutions into four types ac- 
cording to the behaviors of the current densities as follows: 
Type (a) - solutions with < kq < K s , Type (b) - solutions 
with ko ~ kJ, Type (c) - solutions with k\ < Ko < kJ„, 
and Type (d) - solution at k = K m . As we can see in Fig. 1101 
concerning the ratio A4 P /M, the following relations hold: 



> 



M, 

M 
M P 
M 



(Type a) > 
(Type c) > 



M p 
M 
Mp 
M 



(Type b) 
(Typed). 



(60) 



,(9) 



4?r 



(BS 



-Dfl 



<9*(r, 



4-7rr sin 9 dr 



(57) 



Since the solution for the magnetic flux function behaves 
as sin 9, the following quantity becomes a constant and so we 



In FigQj] the distributions of j v normalized by the max- 
imum value of the flux function l^maxj along the normalized 
radius on the equatorial plane are displayed. In this figure the 
interior toroidal current j v , the component of j v due to the k 
term and the component of j v due to the p term are displayed. 
The values of (k , M p /M) in the panels are (a) (3.0, 0.701), 
(b) (~ 4.49, 0.501), (c) (5.0, 0,452), and (d) (5.76, 0.417). 
In Fig[l2] we make that the sign of J to t is always positive. 
In other words, we multiply Jtot/\Jtot\ in order to plot the 
distributions of Jtot- In the left panel of Fig[l2] each line de- 
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Figure 10. Left: Plotted is jr'o against kq for configurations with the parameter satisfying r), ' pop = — 1. We can see that there are two solutions 

without surface currents at ko ~ 5.76 and at kq ~ 9.10 and that there exist two singularities at kq ~ 4.49 and at ko ~ 7.73 within the range 

0. < no < 10.. Right: The energy ratio Ai p /M is plotted against Ko. The solid lines and dotted lines denote the values of kq for the solutions 

without surface currents and the singularities, respectively. The toroidal magnetic field vanishes at kq = 0. 



(a):Ko = 3.0 



(b):K = 4.49 





(c): k = 5.0 



(d):K = 5.76 




Figure 11. The distributions of j<p/c normalized by the maximum strength of Ivfmaxl, i.e. the absolute value of the maximum of the flux function, 
along the normalized radius on the equatorial plane. Each line denotes the value of j v (thick solid line), the component of j v due to the k term (thin 
solid line) and the component of j v due to the p term (thin dotted line). Type (a): In the range of < kq < re* the ip component of the current is 
negative in the whole star and the contributions to the current from the n term and the fi term are also negative. M p I Ai = 0.701 model. Type (b): 
In the range of kq ~ k\ , the contribution to the <p component of the current due to the p term is nearly zero because of the large contribution from 
the K term. M p I M = 0.501 model. Type (c): In the range of kJ; < ko < k^, the if component of the current is positive in most of the stellar 
interior. The contribution to the current due to the k term is positive, whereas the contribution to the current due to the p term is negative. M p I Ai 
= 0.452 model. Type (d): At Ko = k^, the surface current vanishes because jo becomes zero, It is remarkable that not only the component of the 
current due to the k term but also its derivative with respect to the position become at the stellar surface. It corresponds to the null surface current. 
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Figure 12. The total currents normalized l^maxl are plotted against ko- We make Jtot positive in this figure. The left panel: Each line denotes 
Jtot (thick solid line), J sur (thick dashed line), J} J (thin solid line) and J\„J (thin dashed line) respectively. The right panel: Each line denotes 
J tot (solid line) and Jt L (dashed line). 



notes Jtot (thick solid line), J tot (thin solid line), J sur (thick 
dashed line) and J^~ t (thin dashed line). In the right one, we 
decompose the J tot into the total k current and the total fj, 
current. The solid line show the Jf a t and the dashed line show 
the Jt ot respectively. 

It is remarkable that these models have strong toroidal 
magnetic fields contrary to configurations whose magnetic 
fields without surface currents extend to the infinity, i.e. open 
fields, and cannot become large. From the panels in FigfTTI 
the configurations of Types (a) and (b) have the positive sur- 
face current in the <p direction which is opposite to the interior 
negative current density in the <p direction, while the configu- 
rations of Type (c) have the negative ^-surface current while 
the interior (^-currents are positive in almost all of the interior 
region. The model of Type (d) has no surface current. How- 
ever, the if component of the interior currents are negative in 
the finite surface region while the <p component of the interior 
currents in the most inner region is positive. Therefore, in the 
range of < Ko < «J, J tot — Jtot + Jtot an d m tne range of 
kJ < ko < Km, J^ = Jtot and J$~' = J? ot respectively. 

The signs of the k term, the /i term and j v of the solu- 
tions of Type (a) are all negative within the whole region as 
seen from Figllll while the surface current is positive because 
jo > 0. Therefore, the strong toroidal magnetic fields of the 
solutions of Type (a) are sustained by the oppositely flowing 
surface current as has been explained for the reasons why the 
numerically exact open field configurations with surface cur- 
rents can have large toroidal magnetic fields in the previous 
section. For Type (a) solutions, the toroidal magnetic field en- 
ergy increases monotonously as the magnitude of the surface 
current becomes larger. 

The solutions of Type (b) show the extreme features cor- 
responding to the singularity. Since the strength of the surface 
current becomes infinite, the n term becomes larger and the 
contribution from the /j, term becomes nearly zero compared 
with that of the k term. Thus the non-force-free /j, term con- 
tributes essentially nothing to the if component of the current 
density and so those solutions of Type (b) can be considered to 
be almost the same as th ose force-free configurations obtained 
bv lBroderick & Na ravan (2008). 



As seen from these panels, the signs of the /i term and 
the k term are different from each other for configurations 
of Types (c) and (d). In fact, the sign of the surface current 
changes from negative values to positive values between solu- 
tions of Type (a) and those of Type (c). As a result, the interior 
currents can become negative near the stellar surface region. 
Since the surface currents are negative in this parameter range, 
the toroidal magnetic fields can be sustained by both the neg- 
ative surface currents and the negative interior currents near 
the stellar surface region (see Fig|12t. The toroidal magnetic 
field energies become larger in this parameter range than those 
in the parameter range where only oppositely flowing surface 
currents such as solutions of Types (a) and (b) are allowed. 
Moreover, the toroidal magnetic field energy reaches its max- 
imum value for the model (d) which has no surface current. 
These phenomena imply that the effects of the surface cur- 
rents are not so large compared with those due to the interior 
currents near the surface region which flow oppositely to the 
interior currents further inside. 

The model (d) has only the interior negative current re- 
gion without surface currents in addition to the interior pos- 
itive currents in the inner part of the star. The effect of the 
interior negative current is larger than that of configurations of 
Type ( c) (s eeFig|12t. The solutio ns obtained bv lloka & Sasaki! 
(2004) and lYoshida et alj J2012h can be considered to belong 
to the same type as the model (d) except for the compressible 
densities. 

As seen from the left panel of Fig{T2] the total current 
Jtot becomes only slightly larger as ko increases. On the other 
hand, J tot increases rapidly beyond the ko = kJ where j[~ t 
starts decreasing. This means that the negative current region 
can cancel much larger interior bulk positive current than the 
negative surface current. Therefore, as we can see in Fig]12l 
the total k current becomes larger and the total yt current de- 
creases rapidly as ko increases beyond the ko — k\. Since the 
ratio oiMp/M reaches the minimum value for the model (d), 
this kind of configuration without surface currents but with the 
negative interior current region has the strongest toroidal mag- 
netic field energy among all the configurations as far as the 
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functional forms for k and [i are the same as those chosen in 
this paper. 

We have also calculated TV = 1 closed field configu- 
rations. Since we cannot integrate the source term analyti- 
cally for compressible polytropes, we have used the shoot- 
ing method to obtain the eigen solutions for the boundary 
value problems (see lloka & Sasak i 2004). Obtained solutions 
of N = 1 polytropes have the same tendencies as those for 
the TV = solutions. The ratio A4 P /A4 reaches its minimum 
value 0.349 at k^ = 7.42. The toroidal magnetic field en- 
ergy is slightly larger than that of the corresponding TV = 
configuration. 



6 DISCUSSION 

6.1 Open field configurations vs closed field 
configurations 

As we have calculated and seen in this paper, the negative 
(oppositely flowing) surface currents and/or the negative (op- 
positely flowing) interior currents seem to generate strong 
toroidal magnetic fields within the stars. We have obtained the 
configuration having the minimum value of M v j M ~ 0.697 
when jo = 7.5 x 10 -3 for the solutions of open fields by tak- 
ing the surface currents into account. On the other hand, the 
ratio A4 P /A4 reaches 0.349 when all the magnetic field lines 
are closed and confined inside the stellar surface for TV = 1 
polytropes. This value is much smaller than that for open field 
configurations with surface currents. The functional forms for 
the arbitrary functions and the boundary conditions of these 
two models, i.e. the closed filed solution and the open filed 
solution, are different from each other (see Eq|36land EqJ48b. 
but both configurations contain the interior region where the 
interior currents are negative (oppositely flowing) and there 
appear very strong toroidal magnetic fields. 

From the right panels in Fig|7J the signs of the [i term 
and the k term are positive for our open field models with 
negative surface currents. In other words, the models do not 
contain the negative current regions. The open field models 
which we have obtained correspond to models of the Type (a) 
for the closed field configurations. Therefore, the stars could 
have stronger toroidal magnetic field energy if they can con- 
tain the negative interior current regions near the stellar sur- 
faces. However, we could not find the functional form of k 
for which the negative interior current region appears as the 
model of Type (d) in this paper. In any case, the toroidal mag- 
netic field energies of closed field models are larger than those 
of the corresponding open field models with negative surface 
currents. It should be noted that the model of Type (d) sus- 
tains the largest toroidal magnetic field energy among all of 
our solutions obtained in this paper. We can conclude that the 
magnetized equilibria with strong toroidal magnetic field en- 
ergies would be the closed field configurations. 



6.2 Effects of compressibility to toroidal magnetic fields 

We have employed polytropes with TV = 1 and TV = in 
this paper. Since we are mainly interested in the effects of sur- 
face currents, we adopt polytropes as equations of state, al- 
though polytropes are too simple equations of state. The differ- 



ent equations of state result in different density distributions as 
well as different magnetic field structures, because the current 
density formula contains the density depending term (EqjlOll. 
As we have described in Sec l5.2l configurations for polytropes 
with TV = 1 can sustain slightly larger toroidal magnetic fields 
than those with TV = 0, i.e. the minimum value of A4 P /M is 
0.417 for the TV = polytropes and the minimum value of 
M p /M is 0.349 for TV = 1 polytropes as far as the other 
parameters are the same. Therefore, the configurations with 
softer equations of state can sustain stronger toroidal magnetic 
fields for polytropes. 

How about f or more realistic equation s of state discussed 
by other authors ? lKiuchi & Kotakj (200a) calculated twisted- 
torus magnetized equilibrium states using some realistic equa- 
tions of state at zero temperature. Their method is the same as 
our method. Fig. 4 - Fig. 7 in their paper show the density con- 
tours and the magnetic field contours for different equations 
of state. The structures of the poloidal magnetic field lines and 
the regions of the toroidal magnetic fields within the stellar 
surfaces are different among the different equations of state. 
For example, among models of the Shen's equation of state 
(EOS), the position where the toroidal magnetic field attains 
its maximum strength is located near the stellar surface and the 
width of the region where the toroidal magnetic fields appear 
is relatively small, but for models with the FPS's EOS, the po- 
sition of the maximum toroidal magnetic field is shifted to the 
inner part of the star and the size of the toroidal magnetic field 
region is much bigger than that for the Shen's EOS. Although 
they did not calculated the toroidal magnetic field energies and 
the ratio A4 P /A4, they showed the ratio of the local strength 
of the toroidal magnetic field to the poloidal magnetic field (h 
in Table 4). The value of h for FPS's EOS is about as twice as 
that for Shen's EOS. Therefore, although the influence of the 
equation of state might become more important if we would 
consider structures of neutron stars, it would not change the 
basic properties discussed in this paper dramatically, although 
the values and/or the regions for the toroidal magnetic fields 
would surely be somewhat different from those obtained in 
this paper. 



6.3 Stability of configurations with oppositely flowing 
^-currents within and/or on the stellar surfaces 

It is very interesting and important to analyze stability of our 
models for open magnetic fields with surface currents. Since 
some of our solutions satisfy the Braithwaite's stability cri- 
teria, Eq.lQJ, our models could be stable. Although it is very 
difficult to tell the stability for a certain model exactly, we will 
be able to check the stability by several non-exact ways and 
get rough idea about the stability of the configuration. 

First of all, we consider the stabilities of the magnetized 
stars with pure surface currents and with no interior currents 
(see the left panel in Fig.QJ. The stability of the magnetized 
stars with surface currents in the surface region of an infinitely 
thin width could be considered to be essentially the same as 
that of configurations with pure surface currents. If the mag- 
netized stars possess only the surface currents which generate 
the pure dipole magnetic fields outside the stars, their inte- 
rior magnetic fields are uniform along the z axis (see the left 
panel of FiglT}. The magnetic fields of this kind of configura- 
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tion are unstable and decay within a few Alfven ti me, because 
the re is no toroidal magnetic fi elds (iMarkev & Tavlerlll973[) . 
As lFlowers & Rudermanl < l 19771) also explained the ins t ability 
of this kind of configuration and lBraithwaite & Spruit! (2006) 
carried out non-linear evolution of the instability by numeri- 
cal simulations and showed unstable growth of the initial sta- 
tionary states as explained above. Therefore, the fields of the 
magnetized stars with pure surface currents can be considered 
to be unstable. 

By contrast, as for the configurations with surface 
currents, which might lay in, for example, the crusts of 
the neutron stars, the magnetic fields could become sta- 
ble. In such configurations, we can assume that the widths 
of surface current layers are not infinitely thin any more 
and the finite Loren t z forc e acts on the surface currents. 
iFlowers & Rudermanl l l 19771) considered configurations with 
surface current layers as well as with uniform magnetic fields 
and dipole magnetic fields inside and outside of the stellar 
surface, respectively, and found that those configurations with 
current layers might be stable. In realistic situations as neutron 
stars, when the solid crusts of neutron stars form after their 
proto-neutron phase, the crusts could sustain the Lorenz force 
to themselves and they could prevent growth of the instability 
of magnetic fields. For such situations, the magnetic fields can 
survive in much longer time than the Alfven timescale. 

Concerning direct computations of the evolutions 
starting from the perturb ed initial stationary states, 
iBraithwaite & Spruitl ( I2006T) carried out numerical evo- 
lutions of the twisted-torus interior magnetic fields with 
solid crusts. They included surface current layers with finite 
widths as their boundary condition for the magnetar's crust 
and used one of their quasi-stationary twisted configurations 
which they had obtained after long time simulations as initial 
values. Their numerical model is similar to our solution with 
surface currents. The magnetic fields of such stars do not 
decay within the Alfven time scale in their simulations as 
far as the crusts can sustain the Lorentz force. Therefore, 
our twisted- torus models with surface currents would be also 
stable configurations. 

Evolutions and stabilities of configurations for closed 
magnetic fields were argued by iDuez et al.l ( 201 Of). They per- 
formed numerical simulations using IDuez & Mathisl ( 120 101) 
solutions as their initial states. They concluded that models 
with closed fields both with poloidal and toroidal magnetic 
fields do not show any sign of becoming unstable within their 
simulation time if the initial model satisfy the stability criteria 
in Eq. {T). Therefore, the closed magnetic field models which 
are obtained in this paper would be stable configurations. 



6.4 Application to magnetars 

It is important to find out natural mechanisms to generate sur- 
face currents and/or their origins if we apply our models with 
surface currents to real bodies such as magnetized neutron 
stars, especially to magnetars. Magnetars are young neutron 
stars with very strong magnetic fields. The magnetars are con- 
sidered as source objects of special high energetic phenomena 
such as the anomalous X-rays emission and the soft gamma- 
ray emission. Thus those pulsars are called the anomalous X- 
ray pulsars (AXPs) and the soft gamma-ray repeaters (SGRs), 



In particular, their high X-ray luminosities and giant flares 
have been considered to be deeply related to the strong mag - 
netic fields of the stars JThompson & Duncan! 1 19951 120011) . 
The magnetic fields are nearly dipole poloidal fields globally, 
but there would be higher order (such as quadrupole and oc- 
topole) poloidal magnetic fields near the surface or toroidal 
fields winded up by ra pid differential rotation durin g the proto- 
neutron star stages dDuncan & Thompson! I1992D inside the 
star. Before we apply our models with surface currents to the 
magnetars with strong toroidal magnetic fields, we need to 
clarify or at least have some ideas about origins or forma- 
tion mechanisms for the oppositely flowing surface currents 
or the discontinuity of the magnetic fields on the stellar sur- 
faces. Then, what is the origin of the negative surface currents 
or the negative current region? There might be two possibili- 
ties to explain it. One is related to the crusts of neutron stars 
and the other is related to the magnetospheres around neutron 
stars. 

Since the physics of the crusts of the neutron stars is too 
complicated and difficult to deal with, we only assume that the 
crusts consist of highly conductive solid matter. If the crusts 
are highly conductive, the electric currents can exist within 
the crust regions. Then, the crusts can make parallel compo- 
nents of magnetic fields discontinuous near the stellar surfaces 
by the toroidal currents inside the crusts. The magnetic fields 
are frozen to the matter and fixed to the crusts because of 
their high conductivity. On the other hand, the interior mat- 
ter of the magnetars is not solid. Thus the matter inside of the 
crusts can move differently from the crusts and the discontinu- 
ities of the magnetic fields would be born between the crusts 
and the interior regions. The interior fields begin to spread to- 
ward t he stellar surfaces by the so me kind of magnetic dif- 
fusion ( IBraithwaite & Spruitll2006l) the discontinuities would 
be e nlarged by the magne t ic pre ssure. As we have seen be- 
fore, Braithwaite & Spruit (2006) simulated this kind of con- 
figuration and found the growing of the Lorentz stress in the 
crusts. From the direction of the discontinuity, we expect that 
the stress is tensile one globally. If the crusts are cracked by 
the stress, it would result in flares of SGRs. Following this 
scenario, our models with strong toroidal fields as well as sur- 
face currents are considered as stationary states of the crusts 
with strong Lorentz forces before occurrence of giant flares. 
If a part of the crusts is cracking, the magnetic energy and the 
helicity are inje cted from the stars and would produce mag- 
netize d flows JTakahashi et al.ll2009L 1201 ll. iMatsumoto et al.l 
l201ll) . These kinds of magnetized outbursts would be giant 
flares of SGRs. We will consider this process by using our 
models with surface currents in the following. 

At first, the surface currents in the crusts can sustain the 
strong toroidal fields by bending the poloidal magnetic fields 
as shown by the model with jo = 7.5 x 10 -3 in the right 
panel of Fig|9] When the Lorentz force exceeds a certain crit- 
ical value, a part of the crust begins to crack. We can consider 
this phenomena as decreasing the strength of the surface cur- 
rent, because a part of the conductive matter is disturbed by the 
cracking. We assume that a certain cracking reduces the value 
of jo from 7.5 x 10 -3 to 2.5 x 10~ 3 as an extreme example. 
The surface current with jo = 2.5 x 10~ 3 cannot sustain the 
toroidal magnetic fields any more which the surface current 
with jo = 7.5 x 10 -3 has sustained. The toroidal magnetic en- 
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ergy and/or the magnetic helicity would be transferred out into 
the outside of the star in order to relax to the stationary state 
with jo — 2.5 x 10~ 3 (transition from the right panel to the 
left panel in Fig(9). Through this cracking, the dimensionless 
toroidal magnetic field energy M t changes from 2.49 x 10 -2 
to 9.92 x 1CP 3 according to our calculations. Therefore, about 
60 % of the toroidal magnetic field energy would be released 
during the cracking. Although this is an extreme example, it 
is natural that the injection of the magnetic helicity and the 
release of the toroidal field come from the transition of the 
magnetized equilibria by the phenomena such as cracks of the 
crusts which reduce the surface current strength. 

Another possibility is the effect due to the magnetosphere 
which excites oppositely flowing current densities near the 
stellar surfaces. Colaiuda et al. (2008) discussed the impor- 
tance of the magnetosphere as the boundary conditions for 
both the poloidal and toroidal magnetic fields on the stellar 
surfaces. Our present models and many configurations in other 
previous works assume that the outside of the star is vacuum 
where the current density and the toroidal magnetic fields do 
not exit. However, the presence of the magnetosphere changes 
the boundary condition for the magnetic field. This change of 
the boundary condition would significantly influence on the 
magnetic field configurations as we have seen. 

As for rotation powered pulsars, there have been many 
investigations about th eir magnetospheric pheno mena such as 
pulsar winds (see e.g. iGoldreich & Julianl 1 1969ft . Their rapid 
rotations produce enormously large electrical forces and the 
surface charged layers could not stay in their stationary states. 
Charged particles run away from the stellar surfaces and form 
the pulsar magnet ospheres (see, for the recent particle pulsar 
wind simulations, IWada & Shibatall200ll2oTTh . This charged 
particles would produce the strong currents outside the star 
and the twisted magnetosphere would form. 

The rotational speeds and the strengths of the magnetic 
fields of the magnetars are different from those of pulsars, 
but there would be some kinds of magnetospheres around 
the magnetars. The recent X-ray spectral observati ons show 
the p resence of a ma gnetosphere for the m agnetar (IRea et al.l 
120090 . Very recently, IVigano et al.l (1201 lb computed numer- 
ically a force-free twisted magnetar magnetosphere. They 
treated the stationary state of the magnetized star as an in- 
ner boundary condition for the magnetosphere. We can see 
the various magnetospheric structures by applying different in- 
ner boundary conditions in their paper. If the magnetosphere 
forms the oppositely flowing toroidal currents near the stel- 
lar surface, the magnetized star can sustain the strong toroidal 
magnetic field energy inside the star. However, the details of 
the calculation are beyond the scope of the present study. 



7 SUMMARY AND CONCLUSION 

In this paper we have dealt with the effects of surface currents 
upon the interior toroidal magnetic fields. We have shown that 
the oppositely flowing surface currents can sustain the strong 
toroidal magnetic field energy inside the star both for the open 
and closed field configurations. 

In the open field models, we have found that there is an 
upper bound of the total current of the star for a fixed set 
of parameter values. Increasing the maximum strength of the 



toroidal magnetic field decreases the region of the toroidal 
magnetic field due to this upper bound. Therefore this upper 
bound limits the ratio of the poloidal magnetic field energy to 
the total magnetic field energy. To exceed this upper bound, 
magnetized star needs the oppositely flowing surface currents 
to the interior toroidal currents. The interior current can over- 
come the upper bound and the ratio M p /M decreases sig- 
nificantly, because the surface current counteracts the interior 
toroidal current. 

In the closed field models, we have found that a model 
with an oppositely flowing current region but with no surface 
current can sustain the strongest toroidal magnetic field among 
all of our magnetized stationary states. The negative surface 
currents can sustain the strong toroidal magnetic fields in the 
models with closed magnetic fields. However, the strengths of 
the toroidal magnetic fields for models with negative surface 
currents cannot exceed a critical value even if the strength of 
the surface current becomes infinity. In order to overcome the 
critical value, the negative current region is required. Increas- 
ing the size of the negative interior current region decreases 
the negative surface current. As a result, the toroidal magnetic 
fields become the strongest when not only the negative interior 
current region becomes the largest but also the surface current 
disappears. 

It should be also noted that, although we have not im- 
posed a condition d^ /dr\ _ ( a > = 0, the obtained two eigen- 
functions without surface currents fulfill this condition and, 
moreover, that the values of A4 P /A4 are very small. This im- 
plies that by computing a series of eigen states with surface 
currents as well as with oppositely flowing interior currents we 
could have easily reached an eigen state whose eigen function 
behaves very smoothly for which the role of the toroidal mag- 
netic fields becomes very important. Furthermore, it is remark- 
able that these solutions obtained by considering in the wider 
functional space without no restrictions about the slopes of the 
functi ons correspond to tho se solutions obtained by other au- 
thors JDuez & Matiull2O10l) . 

We have applied the models of open magnetic fields with 
surface currents to explain the strong hidden toroidal magnetic 
fields inside the magnetars. We have considered two possibil- 
ities as the origin of their surface currents. One possibility is 
related to the crusts of the magnetars. Since the crusts are made 
of the solid matter, it could make the magnetic fields discon- 
tinuous at the crusts and the surface currents would appear due 
to these discontinuities. The magnetized stars can sustain the 
strong toroidal magnetic field energy by bending the poloidal 
magnetic fields within the crust zones. 

The other possibility for the excitement of oppositely 
flowing currents inside and/or on the stellar surface might be 
related to magnetospheres which produce the oppositely flow- 
ing toroidal currents near and/or on the stellar surface. This 
kind of magnetosphere would also sustain the strong toroidal 
magnetic field energy inside the star. These models might be 
the key to reveal the mechanism of the giant flares of the mag- 
netars. 
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APPENDIX A: PHYSICAL QUANTITIES AND 
ACCURACIES OF NUMERICAL SOLUTIONS 

Al Global physical quantities 

In order to see the global characteristics of magnetized equi- 
libria, we define several global physical quantities as follows: 



W =2 



i d r 



T=^ I p(m) 2 d 3 r, 



U= I pd 3 ' 



F 



x B 



(Al) 



(A2) 



(A3) 



(A4) 



where, W, T, n and F are the gravitational energy, the kinetic 
energy, the volume integral of the pressure and the volume in- 
tegral of the quantity which is an inner product between the 
position vector and the Lorentz force vector, i.e. which is re- 
lated to the work due to the Lorentz force, respectively. We 
also define the magnetic field energies of the star as below: 



Mt = h ' L ' idV ' 



i 

8^- 



M p = ^T ( B r + B e) dV 



M =M P +Mt- 
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Here, M p , Mt and M are the poloidal magnetic field energy, 
the toroidal magnetic field energy and the total magnetic field 
energy, respectively. We use these values to evaluate the ratio 
of \heMp/M. 



A2 Dimensionless physical quantities 

We display dimensionless forms of other physical quantities 
as below: 
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(A8) 
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region, and the other is f — [1, 100], 8 = [0, n], i.e. for the 
vacuum region. We discretize the stellar interior region into 
mesh points with an equal interval dr and the vacuum region 
into mesh points with an equal interval ds where s is defined 
as below: 

(A18) 
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Figure Al. The virial quantity VC is plotted against the number of the 
grid point in the r-direction (N r ) 



(A9) 



tt 



3 A7vGrip n 



V4vrGpn 



K 




-4. 



* 



Jv> 



V4^G/r e 
A, 



\/47rGr 2 p max 



Jv 



(A10) 
(All) 
(A 12) 
(A13) 
(A14) 
(A15) 
(A16) 



where G, p ma x and r e are the gravitational constant, the max- 
imum value of the density and the equatorial radius, respec- 
tively. 



A3 Setting of the mesh points and accuracies of 
numerical solutions 

In order to examine the convergences and the accuracies of our 
numerical solutions, we use the virial relation as follows: 



VC 



W\ 



(A17) 



Since this value must vanish for exact equilibrium solutions, 
we can estimate the global numerical errors by this quantity. 
This value depends on the mesh size, because we use finite dif- 
ference method by discretizing the computational region into 
finite meshes. We use two computational regions because we 
need to integrate Eq. dA6t even at a far distant place from the 
star. One is f — [0, 1], 6 = [0, 7r], i.e. for the stellar interior 



1 +f 



The region of r [1 : 100] corresponds to the region of s[l/2 : 
100/101]. We use the same mesh number in these two regions. 
As for the ^-direction, we discretize it into mesh points with 
an equal interval (d8). We compute some configurations and 
change the number of grid points in the r-direction N r and 
fixed the number of grid points in the ^-direction, N$ — 513. 
Fig lAll shows the convergence of the VC with increasing N r . 
We see almost the 2nd order accuracies from the convergent 
tendency of solutions from Fig. IA1 1 because we use the 2nd 
order integral scheme, i.e. Simpson's scheme. If we use mesh 
numbers, N r — 513 and Ne — 513, the typical VC values 
are smaller than 1.0 x 10~ 4 ~ 10~ 5 . These values are small 
enough to be able to consider the systems are in equilibrium 
states (see also Hachisu 1986). Thus we fixed the mesh num- 
bers N r — 513 and Ne — 513 during all calculations in this 
paper to obtain accurate magnetized configurations. 



APPENDIX B: EXPLICIT FORMS OF K x AND j FOR 
CONFIGURATIONS WITH SURFACE CURRENTS OF 
THE DIPOLE TYPE 

From the Eq.d56t. we obtain the explicit form of K\ as 



K\ = 47rp u p- 
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The coefficient jo for the dipole type distribution of the surface 
current treated in this paper can be expressed as follows: 
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